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Bubbles of inviscid fluid surrounded by a viscous fluid in a Hele-Shaw cell can merge and break- 
off. During the process of break-off, a thinning neck pinches off to a universal self-similar singularity. 
We describe this process and reveal its integrable structure: it is a solution of the dispersionless limit 



CN ' of the AKNS hierarchy. The singular break-off patterns are universal, not sensitive to details of 

the process and can be seen experimentally. We briefly discuss the dispersive regularization of the 
^/J , Hele-Shaw problem and the emergence of the Painleve 11 equation at the break-off. 
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I. INTRODUCTION 

A Hele-Shaw cell is a narrow gap filled with an incompressible viscous liquid. When another incom- 
pressible fluid with low viscosity is injected into the cell, two immiscible liquids form a sharp interface 
which can be manipulated by pumping in or out the fluids. The motion of the interface is called Hele-Shaw 
flow and it belongs to a large class of processes referred to as Laplacian growth (LG) ||l|]. The latter de- 
scribes the evolution of two dimensional domains where the normal velocity of the boundary is proportional 
pj . to the conformal measure — the gradient of the conformal map of the domain taken at the boundary. The 

dynamics of the interface is unstable; An originally smooth interface tends to evolve into a complicated and 
highly curved finger-like shape consisting of a collection of near-singularities |2]. 

Recent progress in the study of Laplacian growth is due to the uncoveringthe integable structure of the 
problem and linking to the universal Whitham hierarchy of soliton theory 
the integrable structure of the problem had appeared earlier in the seminal works ||9 , 

integrable structure brings the Laplacian growth into a wider context of mathematical problems appearing 
in the theory of nonlinear-waves, and recently in random matrix theory and string theory. The relations 
between the integrable hierarchy and Hele-Shaw flow are transparent and geometrically illustrative and, 
therefore, the Hele-Shaw flow may serve as an intuitive geometric interpretation of the formal algebraic 
objects appearing in the theory of non-linear waves and string theory (for a recent development see 113]). 
The relations read: (i) the real section of the spectral curve — the central object in Whitham theory — happens 
to be the interface in the Hele-Shaw flow, and (ii) the cusp-like singularities of the interface are linked to 
shock-waves of dispersive non-linear waves. 
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FIG. 1: A scheme of an experiment where air is drawn from the air bubble surrounded by oil. A simply connected 
bubble pinches-off and eventually breaks-off through a singular neck. After the break-off both bubbles are maintained 
at the same pressure. The singular shape of the neck is universal; It does not depend on the shape of the original 
bubble. 



Soliton theory provides a framework for the study and the classification of singularities emerging in 
Hele-Shaw flows. A universal properties close to a separated singularity is described by a certain reduction 
of a more general hierarchy. For instance, some cusp-like singularities were recently identified with the 
KdV hierarchy laLZD and modified-KdV hierarchy |al, and have been studied in this framework. 

In this paper we study another class of singularities in the Hele-Shaw flow: those arising at the break- 
ojfs — the processes where one bubble with low viscosity breaks into two. We describe the break-off and 
show that it is a solution of dispersionless limit of the AKNS hierarchy (dAKNS). Curiously, the breaking 
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interface has emerged as a target space in the minimal superstring theory 

As a natural extension of the result, we discuss the dispersive regularization of the singularities and 
the emergence of the Painleve II equation. In the next section, we suggest an experimental set-up where 
the surface tension and other suppressing mechanisms of singularities are negligible, so that the universal 
singular behavior at the break-off can be achieved experimentally. 



II. SET-UP AND RESULTS 



A Hele-Shaw cell is a narrow gap between two plates filled with incompressible viscous fluid, say oil. 
Another immiscible and relatively inviscid fluid, say air, occupies part of the cell forming one or several 
bubbles. Most commonly the flow is produced when air is injected into the cell through a single point in the 
cell, pushing oil away to the border of the cell. Here we will study a different process: air is extracted from 
the air bubble — that already exists in the cell — from two fixed points of the cell situating in the bubble as 
depicted in Fig.Q As the air is drawn out, a simply-connected air bubble will break into two. The set-up is 
designed such that the two separate bubbles after the break-off maintain the same pressure. This condition 



is important. 

The Hele-Shaw flow is driven by D' Arcy's law: velocity in the oil domain is proportional to the gradient 
of the pressure: 

v = --—VP. (1) 

The coefficients are the separation of the plates, b and the viscosity of oil, fi. Below we set the pre-factor 
^- to |, and set the flux-rate such that the area of the air bubble is vrt where t is the time of the process 
in the traditional Hele-Shaw flow; t flows backwards in time in our experimental set-up. As the oil is 
incompressible the pressure is a harmonic function in the oil domain. Inside the air bubble the pressure is a 
constant (which may be set to zero) in time and space. If the surface tension is neghgible then the pressure 
is continuous across the boundary. These set up a Dirichlet boundary problem: 

AP = in oil, P = on the interface, P^— log|2;| a.tz—>-oo, (2) 

where we followed the conventional setup of injection, z = x + iy is the complex coordinate. 

A comment about surface tension is in order. In the case of an expanding interface, when the air is 
injected into the cell, almost every initial shape of the interface results in a cusp-like singularity at some 
finite time, if the surface tension and any other damping mechanism is neglected. After this time the 
jroblem defined by Eqs.QandlSlis ill-posed. Introduction of a damping mechanism (such as surface tension 
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1711 ') usually ruins the analytical and integrable structure of the problem. The only exception, where the 
singularities are curbed, is the situation proposed in recent papers f6,r\. 

However, the reverse process, which we consider in this paper, is well defined (section 7 of [6]). Before 
the break-off the interface remains smooth regardless of how thin the neck is. Remarkably the problem 
remains well defined even after the break-off when the tips of the two emergent bubbles are highly curved. 
This is true when the two bubbles are kept at the same pressure, which motivates the design of the cell as 
in Fig. n In this case, the surface tension can be safely neglected by increasing the extraction rate. This is 
perhaps the only case that the singularities can be regularized without resorting to a damping mechanism 
and thus preserving the analytic structure of the zero surface tension problem. 

We will show that the singular behavior of the neck at the break-off falls into universal classes character- 
ized by two even integers {p, q) = {An, 2), n = 1, 2, • • • , and only by a finite number (2n) of parameters. 
They are called deformation parameters. While, an infinite number of remaining parameters characterizing 
the interface become irrelevant in the vicinity of the break-off: the break-off is universal. The deformation 
parameters are the. flows of the integrable hierarchy. 




FIG. 2: Typical break-off of two droplets at successive time slices. 




FIG. 3; Local curves of (4,2)-break-off just before and after the break-off. 



In Cartesian coordinates the shape of the neck reads: 



y' = x^Fp^g 



^lApg 



where we translate the time t such that i = corresponds to the break-off point. In this paper we consider 
only {p, q) = {An, 2) with n being a positive integer. We have obtained A4„ 2 = ^- The function F4n,2, 
a polynomial of degree An at most, is characterized by 2n parameters. Two of these parameters may be 
eliminated by re-scaling the coordinates. For example, in the most generic case of (4, 2)-break-off, shown 
in Fig. |3l the function reads 
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(1 + i?^) ' * > • before break-off 
1 - ^2 ^ t<0: after break-off 

The same curve will appear as the spectral curve of the AKNS hierarchy. 



III. GLOBAL SOLUTIONS 



Before emphasizing the singular behavior at the break-off, we start by a general description of Hele- 
Shaw processes from a global standpoint. We start from the case of an arbitrary genus (the genus is the 
number of bubbles minus one). Then we specify it to the case of two bubbles using elliptic functions. In 
this section we employ a traditional approach of describing the evolution by a set of moving poles of an 
analytical function mappin g the exterior of the air bubbles to a standard domain. This approach has been 



developed in 



1311 . It does not appeal to the notion of integrablilty. 



A. The Schwarz function 

Let us consider a bounded, multiply connected domain D in the complex plane (a Hele-Shaw cell). The 
domain is occupied by an incompressible fluid with negligible viscosity (air bubbles). Unbounded exterior 
of the domain, C \ D, is occupied by an incompressible fluid with high viscosity (oil). Each air bubble 
contains a source through which air can be drawn out of (supplied into) the bubble such that all bubbles are 
kept at an equal pressure. Another neutralizing source is placed at infinity. 

We will assume that the domain D is an algebraic or quadrature domain flSf], or, more generally, an 
Abelian domain 11411 . These domains are defined by the nature of the Cauchy transform. 



Hz) = - f 
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h+(z) at z gD 

(3) 



z + h-{z) at zG C\D 

For quadrature domains the function /i+ is rational: containing a finite number of poles. For Abelian 
domains dh^ is rational and h'^ may have logarithmic singularities. The use of Abelian domains is not very 

o 

restrictive; Any domain with a smooth boundary is known to be a limit of Abelian domains I19r — ^the set of 
Abelian domain is dense. 

The main property of Abelian domains is: the analytical continuations of the function z from different 
connected pieces of the boundary to the exterior C \ D give the same result. The function being obtained 
by this procedure is called a Schwarz function 12( 1 12111 . In terms of the Cauchy transform ^, 



S{z) = h+{z)-h-{z). 

If the boundary is smooth, the Schwarz function can be analytically continued on some strip inside the 
domain D. 

To fix the notations and emphasize the analytical structure of the Schwarz function in C \ D we write 

K L 

k=l ^ ^^' 1=1 
where ~ stands for an asymptotic equality near each of the singularities q^ G C \ D for A: = 1, • • • ,K + L. 

/ifc's are positive integers. 

The evolution of the domain is determined by the holomorphic function, 

4>{z) = ^{x,y) +ip{x,y), 

where ^ is stream function and p is the pressure. In terms of cl){z) and the Schwarz function, D'Arcy's law 
reads, 

dtSiz) = id<Piz). (5) 




FIG. 4: Two distant points A and B along the boundary become close at the break-off. A blow-up of the critical 
break-off region is plotted. Dotted lines are the interfaces when the bubble breaks up. 



An immediate consequence is that the Schwarz function behaves at infinity as, 



/^o 



dtS{z) ~ ^0 "I at z ^ oo, g is constant and dt^i^ = 1. 



(6) 



which follows form the boundary condition dcl){z) -^ —i/z at infinity. Another consequence is: /i+ does 
not depend on time (the proof is in Appendix A): 

dth+{z) = o. 

Therefore the parameters in Q are fixed and can be used as the initial data for the domain. 

For a single bubble, these data are sufficient to determine the evolution. For a domain consisting of 



Q; 



g + l bubbles, additional g parameters are required |5!]; they are, for examples, relative areas of the bubbles, 
pressure differences in the bubbles, or a mixture of the two. Although later we will focus on the equal 
pressure case, here we do not demand this. 



B. Dispersionless string equation 



D'Arcy's law written in ^ acquires another important form if the univalent function cp rather than z is 
chosen as a coordinate: 



dS_d^_dS_dz_ __. 
del) dt dt d(j) 

Here the derivative in time is taken at fixed (p. This form is known as the dispersionless string equation. 



(7) 



A Cartesian form of this equation is convenient to study break-off. Let us choose a coordinate system 
where the origin is at the pinch-off point and the x axis is tangent to the interface at the pinch-off (the dashed 
lines in Fig.|4li. It is convenient to use the following variables. 

iiz)^^{z + S{z)), y{z) ^ ^^{Z - S{Z)). 

Note that x and y are not Cartesian coordinates of the physical plane; They are holomorphic functions 
defined outside the bubble and analytical in some vicinity of the boundary. At the interface, however, x and 
y are equal to its Cartesian coordinates. Using the new variables (with a rescaling of time), Eq.Qis written 
as 

., ,, dydx dydx 

And equivalently, 

dty{x) = -dx4>{x). (9) 

This form of the evolution equation appeared long ago in 12211 . It was recognized in IJ] as the dispersionless 
string equation known in soliton theory. 

C. Hele-Shaw flow as evolving uniformizing map 

It is illustrative to reformulate the Hele-Shaw problem as the evolution of a map z{u) uniformizing the 
exterior domain C \ D. This map maps the exterior domain C \ D consisting oi g + 1 holes, onto a Riemann 
surface U+. 

It is customary and convenient to consider another Riemann surface U_ that is obtained from U+ by an 



anti-holomorphic automorphism: u ^ u* = f{u) satisfying 



S{z{u)) = z{u*) = z{f{u)). (10) 

At the interface, since S{z) = z, we get u = u*: the interface is invariant under the automorphism. The 
invariance allows U+ and U_ to be glued along the interface to a Riemann surface without boundary, U. 
The genus of this surface is g — the number of the bubbles minus 1. This Riemann surface will be invariant 
under the square of the automorphism, which is in this case called an anti-holomorphic involution. The 
Riemann surface obtained by this procedure is called a Schottky double 12311 ; We shall refer to this Riemann 
surface as "the mathematical plane". 
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Summarizing: U is a Riemann surface with an anti-holomorphic involution; there are (7 + 1 closed cycles 
(the interfaces) that are fixed points of the involution; U is divided into U+ (an image of is the oil domain 
C \ D) and its mirror U_, that border each other along the cycles and are permuted by the involution. 

The anti-holomorphic automorphism m — > m* on the mathematical plane induces an automorphism on 
the physical plane: 



z ^f z* = z{u*) = S{z), 

which is called Schwarz reflection. The interface is invariant under the map. 

The Advantage of the above construction is that the analytic functions defined in the oil domain, such 
as S, X, y, and 0, can now be extended to meromorphic functions on the Riemann surface U. From this 
standpoint a Hele-Shaw flow is seen as an evolution of a special Riemann surface. A break-off changes its 
genus by one. To describe the break-off it is sufficient to take two bubbles mapped onto a torus (g = 1). 
Later we will illustrate the above construction on this example . 

We normalize the map such that it maps a point uq to infinity. Being a univalent map it has a simple pole 
at u = liQ E U+. 

z{u) - "° noGU+. (11) 

[u - Uo) 

All the other singularities of z{u) must reside in U_. 

From (ITUl and Q, the uniformizing map (of an Abelian domain) has the following singularity-structure, 
compatible with that of the Schwarz function: 

K L 

fc=i ^^ ^^^ 1=1 

where X^^^^ <^K+k = 0, to ensure the univalency of z{u). The algebraic nature of the domain wan^ants the 
number of singularities to remain the same during the evolution. Again, they are all located in U-. 

In this language the Hele-Shaw flow can be seen as the evolution of a Riemann surface, the domain of 
z{u), which in its turn is governed by the motion of the singularities: a^'s and u^'s. In other words the 
Hele-Shaw evolution is reduced to the pole dynamics of the map z{u). 

Eq. Q21relates the singularities of z{u) in ilTh to the fixed singularities of S{z) in ©. The relations are 
explicitly: 

-dz{ul) 



Qk = z{ul), Hk = Oik 



hk 



J'K) 

where we assume h^ = for k > K + I. Similarly, Q and (fTTT) give 



k = l,2,--- ,K + L, (13) 



qo = z{ul), dtfJ-o = dt ao dz{u^)f'{uo) 



1. (14) 



U_o en 



r+i 




— T 



=T+1 



FIG. 5: Mathematical plane (left) and the physical domain (right). Conformal map, z{u), maps 11+ to the oil-domain 
C \ D. We can see that the grid hnes in U+ are distorted in C \ D by the mapping. The thick lines at Rc(w) = and 
Re(w) — 1 correspond to the interfaces. The red line at the top (bottom) maps to the critical region. 

In the appendix we discuss the derivation of these equations. Together they give a set of algebraic equations 
completely determining the flow in the case of one bubble. 

Additional g constraints will specify either the areas of the bubbles or pressure differences. The area 
of a bubble is the integral of the Schwarz function along the interface: f S{z)dz, where a is the cycle 
encircling the bubble. Pressure-difference between two bubbles, say the ith and jth bubble, follows from 
dSj. It is the real part of the integral of the velocity dtS(z) along the 6-cycle bjj, comiecting ith and jth 
bubble. 

Im0(z)|^. = Re / dzdtS{z) = -dt I (I Siz)dz] . 

In the following section, we choose the pressures equal, as required by the set-up in Fig. ^ We would hke 
to emphasize that the evolution with the same pressures differs drastically from the evolution with different 
pressures, which can also be achieved by controlling the rates of extraction between two bubbles. We don't 
consider such cases here. 

D. Two bubbles 



A Riemann surface describing two bubbles, has the topology of a torus. The torus is modelled as the 
complex plane modulo the periods 2a;i and 2a;2- The anti-holomorphic involution is chosen as f{u) = —u, 
and requires lui to be real and uj2 to be imaginary |14]. Choosing a proper scale we set cui = 1 and 
UJ2 = T, where r is an imaginary number — the modulus of the torus — which changes with time. We take 
the fundamental domain of U as — 1 < Re(u) < 1 and — |r| < Im(u) < \t\ with the opposite edges 
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identified as in Fig. |5| The fixed points of the involution are the lines: Re(n) = and Re(it) = 1. Each 
corresponds to each bound^ of the bubbles. 



idary 
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We will follow Refrs. Il0l.ll4ll to construct a map z{u). For simplicity we discuss only the quadrature 
domains, where there is no log-singularity in z{u). The map can be written in terms of the Weierstrass zeta 
function C{u)- a quasi-periodic function satisfying, 

C{u + 2uJi) = C{u) + 2r]i, r]i=C{uJi), i = 1,2. 

with a simple pole of residue 1 at the origin C{u) ~ u^^. 

The map z{u) reads 

K 

z{u) = aoC{u - uq) + (3 + ^ OjCiu - Uj) (15) 

where uq £ C/+ and n^ G [/_ for j = 1, • • • , K. The condition = Ylj=o ^j ensures the periodicity of the 
map. It reduces the total number of complex variables (excluding the modulus r) to 2K + 2. The constants 
of motion are computed from (fT3l : 

pij = cijdziu*) and qj = z{u*), (j = 1, • • • , K) (16) 

and from (I14t : 

qo = z{uq) and fiQ = 7xodz{uQ) = t + const. (17) 

Hence the total number of equations is 2K + 2 as well. 

The modulus r is chosen to ensure the equal pressure between the two bubbles: = ^dt ^ S{z) dz, 
where b is a cycle traversing the two bubbles. Integrating this equation we obtain another constant of 
motion, 11: 

n = l (f S{z) dz. (18) 

Before a break-off the cycle b does not exist, and at the break-off it contracts to a point. Therefore 11 must 
be zero for a break-off: 11 = 0, which provides the equation to determined the modulus. In terms of our 
conformal map: 

= n = ^ miajQj - UjJIj - UjHj) - ^ ajUk i ^- — + ({uj + Uk)p(uj + Uk) j , (19) 

3=0 j,k=o ^ ^ 

where r is implicit in all the elliptic functions. 

In the conventional set-up where air is injected into bubbles, a smooth merging occurs only for a fine- 
tuned initial condition which satisfies 11 = (IT9l . If 11 7^ 0, one bubble develops a singular cusp before 
the merging. After this moment the problem without surface tension becomes ill-posed, and requires a 
regularization. 
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IV. BREAK-OFF 

To illustrate a break-off we restrict ourselves to a left-right symmetric case. In this case a bubble will be 
broken into two symmetrical bubbles: one is a mirror reflection of the other. The pressure is also symmetric, 
and n is vanishing automatically. 

It is convenient to follow the reversed evolution, that is, to consider a merger rather than a break-off. 
At the merger the period {2uj2) blows up; a rectangular fundamental domain becomes an infinitely long 
channel — a cylinder; the doubly periodic function C,{u) degenerates into a periodic function with period 
2loi = 2. Accordingly, the (imaginary) modulus r gets infinitely large. Therefore, a convenient small 
parameter for an expansion is e*'^^ < 1 (the nome of the Jacobi theta function) which becomes zero when a 
merging occurs. Here we only consider the situation up to the merging point, not after: t < 0. 

The area near the merger corresponds to the domain near the top of the rectangle. In u-coordinate this 
region is pushed away to the infinity as e*'^'^ ^ 0. A new coordinate u = u — t whose origin resides at the 
top of the rectangle is proper. We use the expansion in small e*'^^: 

({u) = C(n + r) = — (1 - 24e2*''^)'u + 27re*''^ sinvrn + Ine'^''^'' sin27ru + • • • , 

up to n-independent constants. The map iT5l expands: 

z{u + t) = C + e*''^ {Ai sin vrn + Bi cos ttu) + e^*''^ {A2 sin 2ttu + B2 cos 2ttu) H , (20) 

where 

K K 

An = 2iTy afcCosnTTUfc, S„ = — 27r> a^sinnvrufc, n = l,2. 
fc=o fc=o 

For a generic choice of the parameters shape of the interface is not physically acceptable: it may contain 

self-intersecting interface. 

A particularly simple configuration occurs by further imposing up-down symmetry. (A break-off with 
no up-down symmetry is shown in Fig.|6l) The full symmetry (up-down and left-right) may be realized by 
setting (i) uq = ^, (ii) the points n^ in U_ are all on the line Re(n) = — ^ symmetrically with respect to the 
real axis, (iii) the residues a^'s are all real, and the residues at poles n^ and Ujt are equal, iv) the merging 
point is at the origin: z{^ + r) = 0. Then we get C = and Ai = B2 = 0. 

Within this configuration let us solve the string equation (0, which is modified into the following form 
using (fTOl : 

dz{u)dtz{—u) — dtz{u)dz{—u) = id4>{u). (21) 
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One use the expansion of z{u) (1^ to obtain the l.h.s. of the above equation d^ : 

(l.h.s.) ~ 27re'''^ -^- — cos vrn, (22) 

where the two remaining constants, A2 and Bi , are reduced to 1 by a rescaling of coordinates. 

To obtain the r.h.s. of d2n i we must find a meromorphic function (p{u) that has vanishing imaginary part 
along the lines Re(n) = and Re(n) = 1 and has the asymptotic behavior of 0(z) = i log{u — uq) ~ 
—i log z near infinity. One can construct such function (l){z) using the Weierstrass sigma function a{u): 

o'iu — Uq) 
cl){u) = i log — — + 2i riiRe{uo) u. 

a{U - Uq) 

The first term renders Im[(/)(n)] constant along each line Re(n) = and Re(n) = 1, then the second term 
shifts it to zero for both lines. In the critical regime one can also expand (/)(u) to produce: 

(r.h.s.) = id(f){u) ~ in e^'^'^ costtu. 

Combined with the l.h.s. (l22l we obtain '^ ^^^^ = const., or equivalently, 



with t = being the merging point {t is negative before merging). 
The evolution of the interface is described at the leading order: 



(x,y)~(y|t| coshvrs, |t|sinh27rs) 



with s € M parametrizing the interface. Alternatively 



This is the simplest and the most generic (4,2)-class of the break-off. 

V. CRITICAL REGIME 

While breaking-off the shape of the interface near the break-off point does not depend on the details of 
the rest of the bubble: locally, the breaking interface is universal. We will call this regime critical. 

In the critical regime, different scales separate. It is convenient to define the scales in the complex plane 
of (j) — the complex counterpart of pressure — which can be put to zero at the the tips. As we did in the 
previous section, it is convenient to think about break-off in a reversed time. Then the process appears as a 
merger. 
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FIG. 6: Up-down symmetry is broken. Thick line is the boundary. Thin lines are equi-pressure lines. The dot is one 
of the stagnation points. There is a merger in the upper part. Another merger in the bottom is imminent, whereupon 
the bubble will form a ring-type domain. 

A critical regime is the image of a small domain containing the point (p = 0. We saw that at </> ~ e*'^'^ we 
see the details of the tips of the merging bubbles and details of the neck after they merge. At 1 » (/> » e*'^'^, 
however, the details of the neck and the tip cannot be seen, and the details of the remote part of the bubble 
can not be seen either. At this scale we see only the asymptote of the neck y^ ~ x^. Finally, (/> = 0(1) 
corresponds to the rest of the body of the bubble. The scale e*'^^ changes with time and eventually disappears 
at the critical point. In that limit, the shape is scale invariant as in Fig. ID Two even integers p and q 
characterize the critical regime. They are sufficient to determine the scaling functions. 

We find solutions only for (p, q) = (4n, 2), where n is an integer. A similar analysis has been carried 
out for cusp-like singularity of Hele-Shaw flow in ||2|]. 

A. (4,2)-break-off 

We start from the most generic (4,2)-class, y oc A±x'^, where the asymptotes may differ at both sides of 
the X axis, A^ / —A-. More degenerate (4n, 2) will be addressed later. 

After break-off we expect two distinct tips, located respectively at x = ui and x = U2. Assuming these 
tips are smooth, we have an expansion: x — Uj oc y^ + 0{y^), which leads to branch cut singularities at ui 
and U2. Together with the (4, 2) asymptote we are led to the following form for the curve: 



y{x) = a{x — b)y/{x — ui){x — U2) + ex + dx + e, (23) 

with b, ui, U2, d, and e vanishing at i = 0. 

The potential <j){x) must have the same branch points as the curve, it must behave as x at large (but 
not too large) x, and is real on the curve (when x is a real coordinate of the curve). These requirements 
determine the potential: 



-\/{x - ui){x - U2) - 7X. (24) 
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The overall normalization can be chosen by a proper rescaling of time t, the coefficient c controls the up- 
down asymmetry of the shape shown in the asymptote y oc (±a + c)x^, the coefficients e and 7 control the 
up-down asymmetry of the velocity. In fact they will be connected by e = jt. 

Putting this ansatz into the evolution equation Q: dty{x) = —dx4>{x), we see that the regular poly- 
nomial part of y and its singular part consisting of the square root are decoupled. The solution for the 
polynomial part set d = 0, and c as a constant, and confirms that e = 7*. 

In order to match the singular parts of dyy{x) and —dx4>{x), it is sufficient to match their large x- 
expansions. y expands as: 

a[x -[h+^^)x + ^ + — + ...J, (25) 

while the expansion of the singular part of —dx4'{x) is: 1 + 0{x~'^). Comparing, we conclude that the 
coefficients of the x^, x and x~'^ terms in d25l are constants and that the constant term is proportional to 
time. The constants must be chosen such that the shape is y = (±a + c)x^ at time t = 0. This is equivalent 

^_^_ / 21/1 

to setting the parameter 11 in ( I19t to zero. We obtain ui = —U2 = V a ^^^ b = 0. Summing up we find 
an evolving curve (we write it in a self-similar form): 



y = x^ [a\ l + —^ + c + -f-^\ . at i<0 (26) 

\ V ax'^ x^ / 

Now we turn our attention to the flow after merging (before the break-off). It is simpler since it doesn't 

involve any singularities. When the branch points meet to become a double point the curve degenerates. 

The potential becomes 0(x) = itx — 72:. A similar procedure yields the following curve as the one before 

the break-off: 

y = x^i±ail + —^\+c + -i-^\. at t>0 (27) 

Equations d^^ and dTTt form our prediction of the self-similar dynamics for the generic (4, 2)-break-off. 

B. {An, 2)-break-offs 

There is a larger class of solutions of the Eq. [Ojfor the same form of the potential (l24l . They belong to 
the (4n, 2) class. Below we only consider the curve with up-down symmetry: symmetric w.r.t. x-axis. In 
this case the curve and the potential read: 



y = P{x)^/{x - ui){x - U2), 4) = -^/{x-ui){x -U2), (28) 

where P{x) is a polynomial of degree 2n — 1. Its zeros are called the double points. Eq.|28l represents a 
degenerate hyperelliptic curve, where all but one branch cut shrink to double points. Double points and the 
branch points ui and U2 move in time as determined by ©. 



15 



m 



A formally identical problem appeared in the study of dispersionless integrable hierarciiies 
Below we adapt some of the results. The idea is based on the following observation: the large x-expansion 
of y in d28l starts with 2n terms of positive powers in x, while the expansion of d(j) starts from the constant 
— 1 and skips the x~^ term. These expansions match only if the coefficients of the positive powers of y do 
not change in time, while the zeroth order term of the expansion of y is equal to t, and the coefficient in 
front of x^^ is a constant. Once these terms are matched the negative powers match automatically. In the 
previous section we saw how it works on the simplest (4, 2) example. 

The functions Hk{x), defined below, form a useful basis of functions with proper analytical behavior: 
these functions have the same branch points as (t>{x), behave as x'^ at large x, and lack all other positive 
powers in large x-expansion: Hk{x) = x^ + 0(1) for fc = 1, 2, • • • . These properties uniquely determine 
them to be given by the formula: 

Hk{x) = I , "^ =1 V(^ - ^Y - ^^ (29) 

_ ~k _ -^fc+i _ Vfc+i 
r x 

where the subscript + denotes the polynomial part of the large x expansion. We set 

V = , r = . (30) 

The two next leading coefficients Rf^ and Vfc are homogeneous polynomials in r and v of degree k which 
may be expressed through Legendre polynomials as 

R^^-^ = {-i)''p''+^coshipPkiismhip) (31) 

y,+i = (-i)'=-V'=+i£^P^(zsinh(^), 

where 

p = r — V , sinh (/?=—, cosh ip = —. 

P P 

A general solution is written as a superposition of H^ with coefficients t/^'s (that will turn out to be 

time-independent) : 

2n 

y{x) = Y,{k + l)tkHk{x). (32) 

fc=i 

The coefficients are the deformation parameters. The differential form of these equation is: 

dtkV = dxHk+i- 
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Matching the leading powers of the expansion ( E^ we get 

dttk = 0, for A; = 1, • • • , 2n 

tr = -Y,{k + l)tkRk+i, (33) 

fc=i 

2n 

fe=l 
where 11 is a constant of integration. The 2n+ 2 conditions d^ determine 2n coefficients of the polynomial 
(l28l and 2 branch points. Thus solving these algebraic equations is equivalent to solving the Hele-Shaw 
problem near the break-off. The last pair of these equations is called the hodograph relation. 

The constant, 11, is the coefficient in front of x~^ term of the expansion of y, and is exactly the integral 
dTSl over the b-cycle of the curve y{x) d28t - the integral along the cycle encircling the two branch points. 
It must vanish for a break-off. 

A particularly simple case occurs at u = when the curve is left-right symmetric. Then 11 = 0, and 
hodograph equation (I33t becomes: 

1 



t = ^(2fc + l)(-l)'=+M '\t,kr'\ 
k=i \ ^ / 

which determines the time dependence of r. 

C. Various break-offs : self-similarity 

Self-similar solutions are of a particular interest. There the evolution of the curve can be viewed as a 
scaling transformation of the coordinate y and x. They occur when all the deformation parameter are set 
to zero except one from even order which we normalize as t2n = 2n+i ■ ^^ ^ consequence of the merging 
condition 11 = we have v = 0, and r remains as the only scale. 



™2n 



y{x) = H2n{x) = ^^ ^5? - r2, (34) 

V va; — r-^ J j^ 

The curve enjoys an explicit formula in terms of the Chebyshev polynomial of the second kind: 

f^r.(coshe) = ^H^y^^. 

/rN2fc^ (2k\ X 

H2k{x) = 2sinh^(-j /^ t^2fc-2j-i(cosh6'). where cosh6' = — 

3=0 \J I 



17 




FIG. 7: (8,2)-break-off at three successive time slices. Only t4 is turned on. The curve is y = x{x'^ — \)\/. 
It shows more singular behavior than the (4,2)-break-off: r ^ \t\^^^- 



The hodograph Eq.|33] degenerates into tr = — i?2n+i|^=o ~ (—1)"'''^ ( | r^"+^ and yields 

n 



r (X \t 



l/2n 



As an example the Fig. illustrates (8, 2)-break-off. 
Another interesting case occurs at U2 ^ cjo while 
the bubbles is so large that one sees only a tip of one bubble. This case has been studied in 



Another interesting case occurs at U2 ^ cjo while ui remains finite. In this case the distance between 



VI. HELE-SHAW FLOWS AND DISPERSIONLESS INTEGRABLE HIERARCHY 



Here we summarize features of the AKNS hierarchy and its relation to break-offs in the Hele-Shaw flow. 
We will show that the flow is identical to the evolution of the spectral curve appearing in the dispersionless 
AKNS (dAKNS) hierarchy. We start by reviewing the AKNS hierarchy 1261 12711 . 



A. AKNS hierarchy 

The AKNS hierarchy is an infinite set of non-linear differential equations written for two functions r 
and V, which depend on time to and an infinite set of deformation parameters (flows) t^'s 1341 . An efficient 
way of defining the hierarchy requires the definitions of a few objects, first. 

The first is the Lax operator L, a 2 x 2 matrix consisting of a differential operator hOt^ (Later we will 
consider the limit, ?i ^ 0, so we keep h in some of the formulae in this section), 13311 

I I iv r 

L = ndt, + u, u = -[ 

\ r —iv 
The second object is the resolvent Q{x), defined as the diagonal component of the Green function: 
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^-fc 






Q{x) = i{tQ\{L — ix^) ^\to), which is formally expanded in the spectral parameter x as 

oo 

Q{x) = "^Qkx' 

k=l 

where the series of 2 x 2 matrices Q^'s are determined from the commutation relation, 

0, 
which, being expanded in powers of x, yield the recurrence relation for Qf^: 



Q{x),L-ix— 



(35) 



2 [^3, Qk+i] = ^OtoQk + [U, Qk]. 



(36) 



The first few terms are given as 
Qi = - 

Q2 - 

Qs -- 



0"! 



0-1 ,a2 , .r 0-3 

-rv r h I , 

2 2 2 2' 



,1 



//nO"1 



,'^2 , .2 '^3 



-(— r + t; r — r") {2r'v + rt>') h zr f- 



where " ' " means the time derivative: hdtQ. For general k, Qk is a polynomial of degree k inr and u, and 
their time derivatives. It follows, from the above definition of Qi, that the equation hdtoL = [Qi — iv^, L] 
holds trivially. 

The dependence on the deformation parameters can be introduced accordingly. Including the time in the 
set of parameters we now define a hierarchy of flows as 13311 



hdt^^L = [Qk+i -ihdt^P—,L], k = 0,1,2,- ■■ . 

where P ^ ji f v dtQ. Eq.|35lensures compatibility of the equations: the flows commute. 
The first two flow equations are given by 



h : ±ihdt,{re^''^) = {re^''^)" -—{re^'f^) 



(37) 



(38) 



t2 ■■ < 



hdt^r = '^r'^r' + Sr'i;^ + ?,rvv' - r'" 



hdt.p = Ir^v + v^- y-^v' - y-^v - v" 

The recursion relations d^ and the flow equations dTTl are more explicitly written in terms of the 
components of Q^, 

Qk = ~^fcy + ^fcy + ^i^-^y • 
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They are: r6fc = — Vf., and 



TZk+i = rVk + vlZk - [^) , (39) 

n+i = ^n + ^n, (40) 



9..r = ^, a.,.= (^). (41) 



and the flow equations d37b : 



It follows that the flows can be written in an integrated form 

= ^(A; + l)tfe7^fc+l (42) 

fc=0 

-n = ^(/c + i)tkVk+i 

fc=0 

with a constant IT. Alternatively, this hodograph-type relations (l42l may be used as a definition of the 
hierarchy. The reader may notice that the structure of these equations is similar to the hodograph equations 
in(|33j. 

There are several known reductions of AKNS hierarchy. 

- NLS: When r is purely imaginary, the first flow equation of the hierarchy is non-hnear Schroedinger 
equation. 

= i?/) + Vj" + —^i^, ip = ±ire"*^, (NLS equation) 

This reduction seems to be irrelevant for the Hele-Shaw problem. Other reductions are relevant. 

- mKdV: When v = and k = ±ir. The even flows constitute mKdV hierarchy. The first equation of 

the hierarchy is the t2-flow: 

3 

= k + -K^K + k'". (mKdV equation) 

where " ' " stands for hdt2 ■ 

- N-AKNS: It occurs if all the deformation parameters vanish except: tjy = 77+1 • ^^^ ^^^^ interesting 
case occurs at A^ = 2. In this case (l42b reads 

^Y+v^r - r" + tor = 
r\ + U = 

which, after ehmination of v yields the (modified) Painleve II equation 

r" = - + ^ + tor. (43) 

2 r-^ 
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B. Linear problem and M operator 



The flow equations appear as a result of compatibility of the spectral problems: 



L^ = ix—^, (44) 



and, 



hdt,^ = (Wfc+i(x) - ihdtj^) ^, (45) 

where Wfc(x) = Qox^ + Qix''-^ + • • • + Qfe. Id 

Another important object of the hierarchy is the M-operator defined as 

M(x)* = a^*. (46) 

Its relation to flows follows from compatibility of the above definition with the Lax equations ( I44l45t 

M{x) = Y,ik + l)tk{nk{x) - ihdt,_,P^). (47) 

The proof is standard in the theory of integrable hierarchy. We sketch it in the Appendix D. 
The relation expressing the compatibihty of d44b and (I46t : 

[M{x),L-ix^] = -i^, (48) 

is sometimes referred to as the string equation. The Af -operator reflects homogeneous properties of the 
wave function ^ with respect to the scaling transformations: tk — > A'^+^t^, and (x, r, v) -^ {j, j, j). 

C. Dispersionless limit and genus-0 Whitham equations 

The Hele-Shaw problem appears as a dispersionless limit of the hierarchy as h ^ 0. 
Let us pursue the limit on a formal level. As is typical for a semiclassical approximations the wave 
function has the form 

^ ~ e"^-^, (49) 

where the "momentum," cj) = 2idtgA, obeys the spectral equation 

det(^^ + U-ix^j =0. 

The latter gives a relation between the "momentum" cp and the spectral parameter 

(/!)^(x) = (x — v)'^ — r^, 
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and determines a spectral surface — a Riemann surface where the wave function is univalent. 

Accordingly, the commutator in (l48l is replaced by a Poisson bracket. If y{x) is an eigenvalue of the op- 
erator M{x), treated as a function on the spectral surface, the pair y{x{(f))) and a;((^) obey the dispersionless 
string equation 

{y,x} = 1. 

An alternative form of this equation follows from P^ . 

dtoV = -dx4> (50) 

where the time derivative is taken at a fixed x. 

The dispersionless limit of Eas. l38l39l43l is easy. One simply drops the time derivative there. In this 
limit the objects of AKNS hierarchy converge to the objects discussed in the previous section fVB] 

T^k — > Rk, 

Vk -^ Vk, 



detlHk-i^^l^Hk. 



Under this identification (H^ become d^^ . and the Hele-Shaw dynamics is seen to be governed by the 
dispersionless AKNS hierarchy. The dispersionless limit can be also understood as genus-0 Whitham equa- 
tions. We briefly discuss it below. For a comprehensive description of Whitham theory see 12911 . 

Integrable equations have special solutions known as finite-gap solutions. They are quasi-periodic in 
times. For these cases, the wave-function ^ as a function of the spectral parameter x is multi-valued. It is 
a univalent on a spectral surface. The case when r and v do not depend on time to is the simplest example. 
The period of this solution is infinite, but the wave function of this solution oscillates with a period ~ /i^^. 
Then the formula ( l49t gives an exact solution. The genus of the spectral surface is zero. 

A class of non-periodic solutions may be obtained from periodic solutions under assumption that periods 
(the moduli) of the solution depend on time but change slowly, such that they do not change much at the 
time of the order to h. Constructions of slowly modulated solutions is the subject of Whitham theory. 

Modulation of genus-0 solutions is simplest. In order to obtain it one simply drops the time derivatives 
in most equations as it has been discussed above. By doing so we obtain the same equations as obtained for 
the break-off in the previous section. 
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D. Hele-Shaw flow and dispersionless limit of integrable hierarchy 

As we have seen, the dispersionless limit of AKNS hierarchy and equations for the Hele-Shaw flow close 
to a break-off/merging are identical. This is not an accident. In |5] an equivalence between Hele-Shaw and 
the universal Whitham hierarchy has been established on general grounds. 

The relation between integrable hierarchies and Hele-Shaw flows can be summarized as follows: in the 
hmit h^ 

- the eigenvalues of the Lax operator become different branches of the potential L -^ (t){x); 

- the operator M becomes the height function M — > y(a;); 

- the real section of the spectral curve y{x) when y = y, x = x is identified with the form of the 
interface; 

- Eq.l48lconverts to D'Arcy's law written in the form ^; 

- The integrable hierarchy describe the response of the interface to a variation of deformation param- 
eters. 

VII. HELE-SHAW PROBLEM AS A SINGULAR LIMIT OF DISPERSIVE NON-LINEAR WAVES 

A relation between the Hele-Shaw problem and the dispersionless limit of integrable equations reveals 
the nature of fingering instabilities and singularities of the Hele-Shaw problem. It is known for a long time 
that the dispersionless limit /i — > of dispersive non-linear waves is singular 13211 . Despite the fact that the 
integrable equations are well defined at all times, the h —^ limit of some — the most relevant — solutions 
cannot be taken after some finite time of evolution. These solutions develop singular features at small scales 
of the order h and have no meaning ath = 0. These singularities appear in non-linear waves as shock waves. 
In Hele-Shaw flows where the surface tension is omitted they appear as cusp-like singularities. 

Some solutions, however, do not lead to a singularity. One of them has been considered in this paper. 
An inevitable break-off of the extracting process in Fig.|2j although singular, is well defined at all time even 
in the limit of zero-surface tension. 

However, an injection process, unless it is fine-tuned to satisfy 11 = 0, leads to cusp-singularities. 

We can summarize this on a simple example of Painleve II Eq.l43] 

^'S^r = !- + \ + tor. (51) 

2 r^ 
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FIG. 8: (4,2)-break-off with nonzero H. One of the tip develops a (2,3)-cusp before merging. Evolution beyond this 
point requires a regularization of the cusp-singularity. 




FIG. 9: Trajectory of the branch points ui and U2 from solving the dispersionless limit of Painleve II equation (15 2t . 
Blue (thick) line is the physical solution and the red (thin) line is the unphysical solution. Being continued through 
the cusp, the solution becomes multi-valued and therefore unphysical. 



As we have mentioned the /i — > limit of this equation, 

^3 Yj2 







+ 



+ tor, 



(52) 



describes the critical regime when all deformation parameters tk vanish for A; > 2. In this limit r and v have 
a simple interpretation: r is a the distance between two tips, u = — ^ is a measure of asymmetry of the 
bubbles (Fig.[Sll. When 11 = the bubbles merge according to the law r = ^— 2to- When 11 / 0, however, 



3n2/3 

' 24/3 



with 



Eq.|52l happens to have four real solutions. The graph of r(to) overhangs at some time t* 
r{t*) = 2^/^n^/^ as in Fig.|9| At this point the tip of one droplet develops a (2,3) cusp singularity y^ ~ x^ 
as shown in Fig. [8l Equation (I52t . as well as the original un-regularized D'Arcy law ^, have no sense 
beyond this point. 

However, the Painleve II Eq. |^ is well defined at any time. The approximation of simply dropping 
the dispersion term h'^dfr is invalid. The reason is simple: an initially smooth function such as r{t) will 
get sharper, such that the dispersion is no longer small. In fact the approximation becomes invalid some 
time before the overhang occurs. This phenomena is known in the theory of non-linear wave as a gradient 
catastrophe. 
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The relation between the integrable equations and the Hele-Shaw problem raises a natural question. 
Can the Hele-Shaw problem be extended beyond ^ ^ limit? Can the parameter h be used to regularize 
singularities in a meaningful physical way? These questions will be the subject of future study. 

VIII. CONCLUSION 

We have analyzed the universal shapes of bubbles in a Hele-Shaw cell at break-off. These are given by 
the Eq. |26l At the pinch-off point the shapes are universal and self-similar: their evolution is reduced to a 
proper re-scaling of the coordinates. 

The universal shapes may be observed in a proposed experiment under the condition that, after break-off, 
the bubbles are all maintained at the same pressure. 

We have extended the relation between Hele-Shaw flow and integrable hierarchies to the case of a break- 
off processes. We have shown that the universal shapes are solutions to the dispersionless limit of the AKNS 
hierarchy and Painleve II equation. 
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APPENDIX A: RICHARDSON'S THEOREM 

Here we prove that h'^{z) does not change in time. The evolution of the domain D ^ D + 5D means an 
evolution of the function as follows: 

6hHz) = -^ f -^ (Al) 

TT JsD Z-W 

where 6D is the area trapped between the two contours dD and d{D + 6D). Its area is given by D'Arcy's 
law 5t \v X dl\ = ^d<p{z) dz, where dl is a tangent line element to the contour dD. This transforms the 
variance of the function (lAll) to an integral holomorphic function 

6h'^(z) = (f> dw4'('w)dw = 

2vr JdD z-w 

which vanishes since the integrand is an analytical function in the exterior domain. 

APPENDIX B: SINGULARITIES OF THE SCHWARZ FUNCTION 



Now we examine the singularities in S{z) (I13ll4t . 
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Let us first pick a pole singularity u^ . In its vicinity S{z) diverges as 



S{z) = z{u*) 



{U* - Uk)h^ ifiu) - f{ul))h^ 



nulf'' {U - Ul)hk 



where we used u = f{u*). In its turn z{u) — z{ul) ~ z'{ul){u — uD as ul locates in U+, where z{u) is 
analytical. Therefore the Schwarz function 



S{z) ~ ak 



f'inl 



1 



{z - zl)h>^ 



l<k< K 



is singular at z^ = z{u'^). Similarly, for the logarithmic singularities yield S{z) r^ at log(z — z^). 
Finally, the singularity (fTTT l at uq determines the behavior of S{z) at z ^ oo 



S{z) = z{u*) r^ z{u^) + dz{u'^){u* - u^) 



z{u*q) + dz{u*o)f (no)(n - uq) ~ z^ + oq dz{u^)f (uq) -. 



(Bl) 



(B2) 



APPENDIX C: BOOST SYMMETRY 



Conventionally, in the AKNS system, the matrix U is defined to be purely off-diagonal 



m, 



2711. It re- 



quires some similarity transformation, called "boost" symmetry 112811 . to obtain our form. The "un-boosted" 
Lax operator is written as 



L = eif^^^Le-i^"^ = hdt„ + 



e*^ 
-iP 



where we have defined l3 = ^ f vdto- 

One may also define the un-boosted Qk as Qk = e'i^^^Qk e~2^'^3 j^jg ^q^ Eq.|^are more natural in 
this frame. 

hdt^L=[Qk+i,Ll fc = 0,l,2,---, 

Also Eq.|45]takes a simpler form in this frame: 

fiatjei^'^^*) = iik+iix) {e-2^''''^) 



sfc-i 



with Hk = Qoi" + Qix^-' + • • • + Qfc. 
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APPENDIX D: 7\f -OPERATOR AND THE FLOWS 

Here we derive the formula 

oo 

fc=i 
from the string Eq. |5H| It is then obvious that, by acting M on the wave-function ^ using Eq. ^5\ one gets 
Eq-El Here, instead, we are going to prove by acting M on the un-boosted wavefunction, ^ = e*'^ 2 ^. 

We first write the M-operator as a Unear sum of all the flows M{x) = YlT=i '^k dt^^^ and prove the 
coefficients Ck to be {k + l)tk. The relation, [M — dx, dt,, — Hk+i] = 0, by acting on ^ = e*'^^~2"^, gives 

the relation 

00 00 

= dxHk+i - dt, ( Y^ cjnj)M + [Hk+i, ( J2 Cfc^j)] 



00 
= = dxTCk+i - Y^ [cj [dt^Hj - [nk+i,'Hj\j + idtkCj)nj 

i=i 

00 

= dxTik+i - Y [^i dj-iUk+i + (9tfcCj)?Yj = d^Tik+i - [diiik+i - (k + l)'Hkj - X](^*fe^i)^i 

= (k + l)nk-Yi9t,c,)n„ 
i=i 

where we have used J2j=i ^j dj-iHk+i = dxHk+i — {k + l)'Hk, following from [R{x), dx — M] = 0. 
This gives Ck = {k + l)tk- 
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